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By Maurice J. Brevoort and Bernard Rashis 
SUMMARY 


TurbuLlent-heat-transfer measurements were obtained through use of a 
three-dimensional axially symmetric nozzle which consists of an inner 
shaped plug and an outer cylindrical sleeve. Measurements were taken 
along the outer sleeve and gave flat-plate results that are free from 
wall interference and corner effects for a Mach number of 5.03 and for 


a Reynolds number range of 5.6 x 10° to 6.5 x 101. Тһе heat-transfer- 
coefficient results are in good agreement with theoretical analyses and 
the recovery-factor results are in good agreement with extrapolations of 
lower Reynolds number date. 


INTRODUCTION 


The design of supersonic aircraft and missiles requires engineering 
information about heat-transfer coefficients and recovery factors that 
extend over a wide range of Reynolds numbers. Eech of the existing experi- 
mental techniques for studying heat-transfer characteristics appears limited 
in various ways. А test setup which utilizes a flat plate for the testing 
surface is subject to corner and edge effects and also shock interference 
if high Reynolds numbers are obtained by extending the plate. 


One existing technique which avoids these effects is that of the 
hollow cylinder mounted centrally in the test section. However, the 
apparatus is subject to the tube choking at low Mach numbers and deviation 
from flat-plate conditions (the boundary-layer thickness approaching a 
sizable fraction of the cylinder radius) when the Reynolds numbers are 
low. In addition, there are differences in the heat-transfer character- 
istics for the inside and outside surfaces of the cylinders. 


The purpose of this paper was to present results obtained by using 
& technique that would avoid the above-mentioned effects and permit sim- 
ple and accurate measurements of local heat-transfer coefficients and 
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recovery factors over an extended range of both Reynolds and Mach numbers. 
The development of this technique for Mach number 5.05 and the results of 
preliminary heat-transfer measurements for the cooling of the surface by 
the airstream are the subject of this paper. 


SYMBOLS 
с specific heat of the sleeve material, Btu/1b-OR 
бо specific heat of the air at constant pressure, Btu/1b-?R 
E acceleration due to gravity, ft/sec? 
n heat-transfer coefficient, Btu/ft*-see-OR 
k | heat conductivity,  Btu/ft-sec-OR 
m mass of the sleeve material,  pounds/ft^ 
Nu Nusselt number, hx/k 
Pr Prandtl number, uCpg/k - 
R Reynolds number, pVx/p 
Teg T T 
ШИ recovery factor, ---------- 
T tag - 1. 
86 | Stanton number, Nu/RPr 
Teq effective stream air temperature at the wall, some tempera- 
ture which gives a thermal potential which is independent 
of conductance h, 
: О 
Т stag stagnation temperature, К 
Tw inside surface temperature of the sleeve, ?R 
T free-stream temperature, ?R 


V free-stream velocity, ft/sec 
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х longitudinal distance along sleeve, ft 

р free-stream density of the air, slugs/ft? 

u dynamic viscosity coefficient, lb-sec/ftf 

dT /àt time derivative of average wall temperature for а particular 


x-station, ОВ/вес 
APPARATUS AND METHOD 


The apparatus consisted of a three-dimensional axially symmetric 
nozzle which was directly connected to the settling chamber of a cold-air 
blow-down jet. The nozzle consists of a shaped plug and an outer sleeve 
ОҒ 8-inch extra-heavy seamless carbon-steel pipe which was machined inside 
and outside to a wall thickness of 0.388 inch. A detailed drawing of the 
apparatus is shown in figure 1. The ordinates of the plug are given in 
table I. The plug was designed by using three-dimensional characteristics 
(see ref. 1). 


Details of the installation of the thermocouples and static-pressure 
orifices are shown in figures 2 апа 5. The thermocouples are located 
0.060 inch from the inner surface of the sleeve. The wires are number- 
30-gage copper-constantan (0.010 inch in diameter). For these wires, 
conduction along the wire length is negligible. As indicated in figure 2, 
the thermocouples are in intimate contact with the steel wall so that 
thermal resistance at this junction is negligible. 


Figure 1 shows the Location of the static-pressure orifices. The 
Mach number distribution was measured both lengthwise and around the 
sleeve. The results are shown in figure 4. 


The temperature-recording equipment consisted of three synchronized 
high-speed Brown recorders, having a total of 36 switches, two of which 
were connected to the settling-chamber thermocouples, the rest being con- 
nected to the sleeve thermocouples. The stagnation-temperature variation 
with time for a settling-chamber pressure of 220 pounds per square inch 
gage is shown in figure 5. Figure 6 shows the wall-temperature variation 
with time for stations 12 and 18 for a settling-chamber pressure of 
220 pounds per square inch gege. Test runs were made for settling-chamber 
pressures of 92, 142, 220, and 290 pounds per square inch gage. Excluding 
the first 20 seconds, the pressures were maintained constant for each run. 
The recorders were calibrated immediately before and after each run and 
are accurate to t1? Е. 
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In figure Т there are plotted for various times during the run the 
curves of wall temperature against longitudinal distance along the cylin- 
der. These curves were used to determine the rate of change of the longi- 


tudinal conduction қат along the cylinder. Test results were taken 


2 
ах 
Ёс: 


only for the length of the cylinder for which 5 was Zero. 
ах 


REDUCTION OF DATA 


Short test runs in which wall-temperature—time histories are 
obtained give the essential data from which heat-transfer coefficients 
and recovery factors may be determined. The reduction of these data 
requires a method in which the stagnation temperature and, accordingly, 
the equilibrium temperature may vary in an arbitrary manner. 


The Stanton number is calculated from 


and the heat-transfer coefficient is calculated from 


h = me ot fat (2) 
Ta 7 Теа 


Ву definition, the recovery factor is 


Цэн = — (5) 
stag | 


The method of reducing the dat& is simply to select а recovery Ғас- 
tor апа from equation (5) obtain Teg: Substitute this Те, into equa- 
tion (2) and obtain values of h for different heat flow rates. These 
values of h are then substituted into equation (1). The true values 


of Teg? Пу» and St are obtained when the St is constant with time 


(different heat flow rates). Figure 8 shows the curves used in evaluating 
the Stanton number and the recovery factor for station 12 for a settling- 
chamber pressure of 220 pounds per square inch gage. 
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The values of specific heat and specific gravity of the sleeve 
material were taken from reference 2. Figure 9 gives the variation of 
the specific heat of the sleeve material with temperature. The values 
of viscosity and Prandtl number (0.71) for air were also taken from 
reference 2 and were based upon the inside surface temperature of the 
sleeve. The value of TT, chosen was at 120 seconds after starting. 


(This arbitrary choice of temperature can be made because the viscosity 
coefficient is practically a constant factor of the heat-conduction 
coefficient.) 


ACCURACY 


The basic data consist of temperature-time histories at each thermo- 
couple location. The recorders were calibrated for synchronization imme- 
diately before and after each run and are accurate to +10 F. The thermo- 
couple wires are number-50-gage (0.010-inch diameter) copper-constantan 
wires. Conduction along the length for these wires is negligible. The 
wires are insulated except at the junction where intimate contact with 
the sleeve keeps thermal resistance to a minimum. 


The sleeve is enclosed in a vacuum jacket to avoid free convection 
and, for the temperature range considered, radiation effects range up to 
2 percent. The product of the mass per unit area and the specific heat 
of the sleeve material as used in equation (2) is accurate to #5 percent. 


The value of GT/dt used in equation (2) is that of the actual 
thermocouple reading, whereas it should be the value associated with the 
average temperature in the thickness of the sleeve at each point. During 
the test runs, the inner-surface temperatures are lower than the average 
temperatures by approximately 1° to 1.5° at the beginning of the runs 
(40 seconds after starting) and become the same at the end of the runs 
(200 seconds). This effect results in the inner-surface temperature-time 
derivative being less than the average temperature-time derivative by 
approximately 1 percent. The computed slopes when integrated reproduce 
the original curves to within 30.59 F. 


The analysis of reference 3 was used to determine the effect of the 
sleeve diameter. This effect was computed to be less than 0.01 percent. 


RESULTS AND DISCUSSION 


Figure 7 shows the variation of wall temperature with longitudinal 
distance along the cylinder for a settling-chamber pressure of 220 pounds 
per square inch gage. Over the test range, the wall temperatures are con- 
stant with x. In the relation 
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me = 
RCM. 
Ta а Теа 


m and с are constant and dT/dt is constant because T, із constant. 
Therefore, if there is to be а variation of h with Reynolds number or 
X, Teg must vary with x. The value of Teg obtained for the run with 


а chamber pressure of 220 pounds per square inch gage, evaluated at 
120 seconds after starting, actually decreases approximately 29. 


Figure 10 shows the variation of local Nusselt numbers with local 
Reynolds numbers. The х used in evaluating these numbers was based 
upon an x that was considered zero at the minimum. A single line 
faired through all the data points has a slope of approximately eight- 
tenths but, if lines are faired through the points for the individual 
runs, the slopes are slightly higher. It is reasonable to assume, how- 
ever, that the higher stagnation pressures would produce earlier transi- 
tion and hence the zero value of x would move upstream with increasing 
stagnation pressure. Hence, if it is assumed that the individual runs 
should merge into a continuous line (rather than a series of steps), x 
may be adjusted for x = О locations which would make all the runs coin- 
cident. This adjustment has been done in figure 11 by using values of 
x equal to zero at 2.00, 4.25, 7.25, and 10.00 inches ahead of the mini- 
mum for the runs with chamber pressures of 92, 142, 220, and 290 pounds 
per square inch gage, respectively. This systematic adjustment with 
pressure of the x = О location does not affect the value of a faired 
curve through all the data (figs. 10 and 11). 


The Nusselt numbers were found to vary from 4,125 to 32,500 for the 
Reynolds number range of 5.6 х 106 to 6.5 x 101 (based on adjusted x = 0 
locations). For comparison, the curves based upon the analyses of Van 
Driest (ref. +) and Eber (ref. 5) are shown. Іп the Van Driest analysis, 
T,/I, was considered to be 5.0 (the average value of the test results). 


For comparison with these references, the data were computed by using the 
free-stream temperature to determine the density and the velocity. The 
wall temperature was used to determine the viscosity and Prandtl number. 
The two experimental values from V-2 data (ref. 6) are also determined 
in this manner. 


The results of this paper are approximately 10 percent higher than 
the Van Driest analysis and approximately 27 percent higher than the 
Eber analysis. The V-2 experimental results coincide with the results 
of this paper. 


Figure 12 shows the variation of local recovery factors with local 


Reynolds numbers. The variation is from 0.887 at 5.6 x 10° to 0.868 at 
6.5 х 101. The results are compared with an extrapolated curve which 


g 


NACA ТМ 3303 T 


represents the fairing of the data of reference T. The data of refer- 
ence Т are for Reynolds numbers below 6.75 x 106, Тһе results of this 
report and those of reference 7 agree to within 0.5 percent. Also 
included for comparison are the curves for recovery factor equal to 


prl/3 and prl/2, The wall temperature was used to determine the 
Prandtl number. 


CONCLUDING REMARKS 


It can be concluded that the method of testing and reducing the 
data gives results that are in good agreement with established theoretical 
and experimental work. The reduction of the heat-transfer data requires 
adjusting the x = О locations (x is the longitudinal distance along 
sleeve) so that they are consistent when the Reynolds numbers are varied 
either by changing the pressure or the "x value" location. The varia- 
tion of Nusselts number with Reynolds number before and after adjustment 
of the х = О locations are essentially the same. 


The recovery-factor results are in good agreement with the extra- 
polated results of lower Reynolds number data. 


Langley Aeronautical Laboratory, 
National Advisory Committee for Aeronautics, 
Langley Field, Va., June 24, 1954. 
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Figure 1.- Test arrangement. 
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Each wire wrapped and covered 
Both wires form one lead from 
each thermocounle. 


Insulation (ceramic) 


Sleeve (0.388 wall thick) 


| х Y (steel) 


Inner surface 
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ido 32 6) 4 Tips of wires silver- 
(steel) soldered together 


Figure 2.- Thermocouple installation. 
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Figure 5.- Pressure tubes installation. 
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Figure 4.- Mach number distribution for settling chamber pressure of 
220 pounds per square inch gage. 
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chember pressure of 220 pounds per square inch gage. 


Figure 5.- Stagnation temperature variation with time for a settling 
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Figure 6.- Wall temperature variation with time for a settling chamber 
pressure of 220 pounds per square inch gage. 
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Figure 7.- Wall temperature variation with longitudinal distance for a 
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Figure 8.- Stanton number as a function of time and recovery factors for 


station 12 for a settling-chamber pressure of 220 pounds per square 


inch gage. 
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Figure 9.- Variation of specific heat of sleeve material with temperature. 


Nu 


Nusselt number, 


3.0 |x 205 


Settling chamber 
pressure, psig 


О 92 
D 142 


© 220 


2.0 


COSE NL VOVN 


ЯЛ .6 2: 1.0 2.0 3.0 4.0 5.0 x 107 
Reynolds number, R 


Figure 10.- Variation of local Nusselt number with local Reynolds number 
for x=0 at station 0. 
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Figure 11.- Variation of local Nusselt number with local Reynolds number 
for corrected x= О locations. Viscosity and Prandtl number determined 


for wall temperatures. 
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Figure 12.- Variation of local recovery factor with local Reynolds 
number for corrected .x = О locations. Viscosity determined 
for wall temperatures. 
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